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ABSTRACT 

m ■ 

fSj ' We investigate the photometric properties of 456 bright galaxies using imaging data 

recorded during the commissioning phase of the Sloan Digital Sky Survey (SDSS). Mor- 
phological classification is carried out by correlating results of several human classifiers. 
Our purpose is to examine the statistical properties of color indices, scale lengths, and 
concentration indices as functions of morphology for the SDSS photometric system. We 
i/-) ■ find that u' — g', g' — r', and r' — i' colors of SDSS galaxies match well with those expected 

from the synthetic calculation of spectroscopic energy distribution of template galaxies 
and with those transformed from UBV Rc^c color data of nearby galaxies. The agree- 
ment is somewhat poor, however, for i' — z' color band with a discrepancy of 0.1 — 0.2 
mag. With the aid of the relation between surface brightness and radius obtained by 
Kent (1985), we estimate the averages of the effective radius of early type galaxies and 
the scale length of exponential disks both to be 2.6 kpc for L* galaxies. We find that 
the half light radius of galaxies depends slightly on the color bands, consistent with the 
expected distribution of star-forming regions for late type galaxies and with the known 
color gradient for early type galaxies. We also show that the (inverse) concentration 
index, defined by the ratio of the half light Petrosian radius to the 90% light Petrosian 
radius, correlates tightly with the morphological type; this index allows us to classify 
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galaxies into early (E/SO) and late (spiral and irregular) types, allowing for a 15-20% 
contamination from the opposite class compared with eye-classified morphology. 

Subject headings: galaxies:fundamental parameters — galaxies:photometry 



2 



1. Introduction 

The Sloan Digital Sky Survey (SDSS; York et al 
2000) has started producing both photometric and 
spectroscopic data. It is expected that the SDSS will 
produce a very large and homogeneous database for 
research on the nature of galaxies. Photometric ob- 
servations are made with the new five color bands (v! , 
g' , r' , i', and z') that divide the entire range from the 
atmospheric ultraviolet cutoff to the sensitivity limit 
of silicon CCD into five, essentially non-overlapping, 
passbands (Fukugita et al. 1996). These passbands 
are chosen on the basis of their astrophysical merits, 
but this specific choice makes it difficult to compare 
or combine the SDSS data with those obtained with 
more conventional photometric systems. 

In this paper, we present results from a brief anal- 
ysis for properties of galaxies using a small set of 
data for bright galaxies in the new photometric sys- 
tem. The present results have already been used 
to construct empirical galaxy models for simulations 
that have extensively been used to tune the SDSS 
project software (Lupton et al. unpublished) and to 
aid galaxy science (Yasuda et al. 2000). 

To study properties of galaxies, it is essential to 
classify the galaxies into morphological types, since 
different morphological types exhibit distinctly dif- 
ferent astrophysical properties, reflecting the differ- 
ent histories of the formation and evolution of galax- 
ies. A number of attempts have been proposed for 
automated classification of morphologies (e.g., Doi, 
Fukugita, & Okamura 1993; Lahav et al. 1996; Abra- 
ham et al. 1994, 1996), but visual classification still 
serves as the most reliable method when we adopt 
the Hubble classification (Sandage 1961) for galaxies 
with large apparent sizes. This visual inspection pro- 
cedure, which is of course very labor-intensive, limits 
the size of our sample. While we are conducting an 
attempt to produce a large-scale eye-classified sam- 
ple, we present here an analysis we have carried out 
as our initial study. 

In this paper, we focus our considerations on three 
statistical quantities: color, effective size, and concen- 
tration parameter. Galaxy colors are an important 
quantity that characterizes stellar contents of galax- 
ies. The colors at zero redshift are also used as a 
fiducial zero point to study statistical properties of 
faint galaxies and their evolution. An important test 
is to confirm whether the colors obtained with SDSS 
photometry are consistent with those expected from 



broad band photometry using the conventional pass- 
bands, B, V, R, and /, and with those calculated from 
template spectroscopic energy distributions of galax- 
ies. This is an initial step in establishing a photomet- 
ric system where the spectrophotometrically synthe- 
sized flux matches closely with the broad band flux, 
which is one of the goals of SDSS photometry. 

The second statistical quantity is the distribution 
of effective sizes of galaxies. The de Vaucouleurs pro- 
file (de Vaucouleurs 1948) is usually characterized by 
the half light radius r e and the exponential disk by 
the scale length h. These parameters set the funda- 
mental length scales of galaxies. How these quanti- 
ties scale with magnitude is indispensable information 
for construction of simulations of galaxies, in partic- 
ular the study of the detectability of galaxies at faint 
magnitudes and the performance of star galaxy clas- 
sification algorithms, both of which depend on the 
apparent sizes of galaxies. It is also interesting to 
ask whether the effective size varies across the differ- 
ent color bands, which provides information needed to 
understand the statistical distribution of stellar pop- 
ulations in galaxies for a global sample. 

The third quantity we study is the concentration 
index of the light distribution of each galaxy. This pa- 
rameter is known to correlate with the morphological 
type (Morgan 1958); earlier type galaxies have light 
profiles more concentrated toward the center. There 
is some evidence that this parameter can be used 
for automated morphological classification (Doi et al. 
1993; Abraham et al. 1994). Specifically, we study 
the performance of automated classification schemes 
when the concentration parameter is used to classify 
galaxies into early and late types. 

We describe briefly the data and the photometric 
catalog in §2. In §3 our morphological classification 
and the resulting catalog are described. Color indices 
are studied in §4. In §5 the distribution of the scale 
lengths of galaxies is discussed, and in §6 the con- 
centration index is studied. Conclusions are given in 
§7. 

2. Observation and photometric catalog 

The SDSS (York et al. 2000) is carried out us- 
ing a wide-field 2.5m telescope, a large format imag- 
ing camera (Gunn et al. 1998), two fiber-fed dou- 
ble spectrographs, and an auxiliary 0.5m telescope 
for photometric calibration. The sky is imaged with 
thirty photometric CCDs, arranged in six columns of 
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five rows, each row corresponding to a different pass- 
band. The data are taken in time-delay-and-integrate 
mode at the sidereal rate along great circles on the 
sky, yielding a strip consisting of six very long and 
narrow scanlines, each 13'. 5 wide. The effective ex- 
posure time is 54.1 sec. The scanlines of a given strip 
do not overlap, but observing a second strip offset 
from the first strip by about 12'. 8 gives a filled stripe 
2.°5 wide. The pixels are 0". 396 square on a side. The 
0.5m Photometric Telescope observes standard stars 
to determine the atmospheric extinction on each night 
and ties the standard stars to objects observed with 
the 2.5m survey telescope. SDSS imaging data taken 
in five passbands, vf, g', r', i' and z' (Fukugita et al. 
1996), are processed with the photometric pipeline 
(hereafter Photo; Lupton et al. unpublished) specifi- 
cally written for SDSS imaging data processing. 

The galaxies used in this paper are taken from ob- 
servations of the Northern Equatorial stripe on 19 
March, 1999 (SDSS Run 745), supplemented with 
observations on 20 March, 1999 (Run 752) and 21 
March, 1999 (Run 756). These imaging data were 
taken before the commissioning of the current Pho- 
tometric Telescope. We calibrated the data by ob- 
serving secondary patches in the survey area using a 
(now decommissioned) 61 cm telescope at the obser- 
vatory site, and by observing primary standard stars 
using US Naval Observatory's 40-inch telescope, with 
filter and CCD characteristics nominally identical to 
those at the SDSS Photometric Telescope. Since the 
transformation from the primary standard stars to 
the objects observed with the SDSS 2.5m telescope 
has not yet been fully defined, we expect photometric 
errors of « ±0.05 mag in each band with respect to 
the proper SDSS photometric system. Thus, in this 
paper, we will denote our magnitudes as u* , g* , r* , 
i* , and z* , to emphasize the preliminary nature of our 
photometric calibration, rather than the notation u', 
g' , r' , i', and z' that will be used for the final SDSS 
photometric system. However, we will use the latter 
notation to refer to the SDSS photometric passbands 
themselves. The magnitudes we use are based on the 
ABgs system (Fukugita et al. 1996). 

The galaxies we study in this paper are taken from 
Run 745, which ranges from RA=10 h 41 m to 16 h 41 m 
covering effectively 127 square degrees in total. We 
take all 217 galaxies brighter than mp{g*) = 16 mag 
contained in Run 745 for morphological classification 
by eye (here mp{g*) stands for the Petrosian magni- 
tude in the g 1 band as defined below). Additionally, 



we take all 160 galaxies with 16 < m p (g*) < 17 in 
the RA range 10 h 41 m to ll h 41 m , and 107 galaxies 
with 17 < m p (g*) < 18 in the RA range 10 h 41 m to 
10 h 44.5 m . This makes the size of our morphologically 
classified sample to be 484 galaxies in total. Our sam- 
ple is not homogeneous with respect to brightness, 
but otherwise we do not expect any other selection 
biases. Those galaxies which fall on the edge of the 
fields are not taken in our sample. The g' band is used 
in order to facilitate comparison of our morphological 
classification with the work done in the past, which 
has been almost exclusively performed using the blue 
band image. 

After we carried out morphological classification 
using Run 745, an improved photometric data pro- 
cessing was done for Runs 752 and 756, the regions 
of which overlap with that of Run 745. Therefore, we 
identified galaxies in our sample with those imaged by 
Runs 752 and 756 and developed a photometric cat- 
alog from these two runs0. Unfortunately, the entire 
region of Run 745 is not covered with Runs 752 and 
756, and this leads us to drop 28 galaxies. Hence we 
work with 456 galaxies in the analysis given in this 
paper. Details of photometry of galaxies in Runs 752 
and 756 are discussed in Yasuda et al. (2000). 

While many photometric parameters are measured 
in Photo, we are concerned in this paper only with 
three classes of parameters: 

(1) Petrosian radius in the five bands: rp(u'), rp(g'), 
rp(r'), r P (i% r P {z') 

The Petrosian radius rp is defined as the parameter 
that satisfies 



' 2nj^I(r)rdr/(nr%Y 

In practice, this implicit equation for rp is replaced 
with 

2 *fo£ P I{r)rdr/[*{{l.2hr P y - (0.8r P ) 2 )] 

Tj = '■ — 5 

2tt J I(r)rdr / {irr 2 p ) 

(2) 

where I(r) is the surface brightness profile of the ob- 
ject. We adopt r\ = 0.2. We measure the Petrosian 
radius in each band, but we use only the radius mea- 
sured in the r' band to calculate the Petrosian flux 



We used Photo version 5_0_3 of late 1999. The most important 
improvement that concerns us here is in the accuracy of fluxes 
of deblended objects. 
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(or Petrosian magnitude) for all color bands, so that 
the aperture for the Petrosian flux is common to the 
five passbands. 

(2) Petrosian magnitude in the five passbands: m P (u'), 
mp(g'), mp(r'), m P (i'), m P (z') 
The Petrosian flux is defined as 



F P = 2ir 



pkrp 

/ I(r)rdr, 
Jo 



(3) 



where k is set equal to 2 and r P is the r' Petrosian ra- 
dius. For further detailed discussion of the Petrosian 
flux, see Yasuda et al. (2000). 

(3) Petrosian half light radius and Petrosian 90% light 
radius: r 50 (Xi), rg (A;) 

r 50 and r 90 are defined, respectively, so that the 
fluxes over the aperture with these radii are 50% 
and 90% of the Petrosian flux in the passband A, = 
u' , g', r' , i', and z' . 

To give an idea about the basic properties of our 
sample, we show the distribution of (/-band Petrosian 
magnitudes in Figure 1, and those of rp(r'), r5o(r'), 
and r9o(r') in Figure 2. Magnitudes are corrected for 
Galactic extinction using the dust distribution esti- 
mated by Schlegel, Finkbeiner, & Davis (1998) with 
the ratio of total to selective extinction R g / — 3.79 
and R r i = 2.75. When necessary we use R u > —5.2, 
Ri> = 2.09, and R z i = 1.48 for the other color bands. 



3. Morphological classification 

The galaxies in our sample are classified into 7 mor- 
phological classes, T — (corresponding to E in the 
Hubble type, or -5 in RC3 type), 1 (SO, T RC3 = -3 
to -1), 2 (Sa, T RC3 = 1), 3 (Sb, T RC3 - 3), 4 
(Sc, T RC 3 = 5), 5 (Sdm, T RC 3 = 8), and 6 (Im, 
7rc3 = 10), where T R c 3 refers to the type index 
defined in the Third Reference Catalogue of Bright 
Galaxies (de Vaucouleurs et al. 1991, hereafter RC3). 
This classification is coarser than the standard RC3 
types, but we consider that our classification is suffi- 
cient for most purposes for galaxy science. Four of the 
present authors (MH, TI, SO, and MS) independently 
classified all galaxies by eye, by comparing the im- 
age of each galaxy displayed on the SAOimage viewer 
with template galaxies of Frei & Gunn (1994) which 
are given morphological types by RC3. The templates 
of Frei & Gunn cover from E to Im. When one can- 
not assign a morphology to a galaxy, an index of —1 
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Fig. 1 . — Distribution of Petrosian magnitudes in the 
g' band for 456 galaxies analyzed in this paper. Mag- 
nitudes are corrected for Galactic extinction. 



is given. For each galaxy, we adopt a median of four 
classifiers as the final value of T (a half integer means 
that the median falls in the middle). In finding a me- 
dian we ignore T = — 1 when one or two classifiers 
give this index (no galaxy was given — 1 by more than 
two classifiers). In this way some galaxies are given 
half integer values of T. The catalog will be published 
elsewhere (Fukugita et al., unpublished). 

We have identified the 54 galaxies listed in RC3 
and in Run 745. Figure 3 compares our morphol- 
ogy T(ours) with that given in RC3, T R C3, for the 
54 galaxies. The size of the symbol is set so that its 
area is proportional to the number of galaxies in the 
grid. The curve connects points where our morphol- 
ogy matches with RC3 morphology. A tight corre- 
lation is seen between the two classifications. For a 
given T(ours) the average discrepancy in the two clas- 
sifications (in units of T R C3) is (|AT R C3|) = 1.6 with 
the rms scatter of 1.8. Nairn et al. (1995) compared 
experts' classifications and concluded that the rms 
scatter among classifiers is 1.8. Somewhat larger dis- 
crepancy between our classification and that given in 
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Fig. 2. — Distributions of Petrosian radii, Petrosian 
half light radii, and Petrosian 90% light radii of the 
456 galaxies measured for r'-band images. 

RC3 may be expected, owing to the fact that we have 
estimated morphology from CCD images for which 
the contrast and grey scale are not fixed, whereas 
RC3 classification is based on photographic materi- 
als. On the other hand, Nairn et al. used the identical 
prints for test classifications. In view of these consid- 
erations we take our classification to be acceptable. 
In specific details, our classification tends to classify 
some SO galaxies as E, and classify spiral galaxies into 
somewhat later types. 

Figure 4 shows the morphological type distribu- 
tion of galaxies in our sample. The relative fractions 
of types given in the literature (Fukugita, Hogan, & 
Peebles 1998; Loveday 1996) obtained in the B or bj 
band are also plotted for comparison (normalized to 
the total number of galaxies in our sample). Love- 
day does not classify spiral galaxies into subclasses; 
so his numbers are distributed evenly over classes 2, 
3, and 4. In our sample the number of galaxies with 
T = (i.e., E) is significantly higher than that of 
T = 1 (SO) galaxies, in contrast to the other results, 



Fig. 3. — Morphology determined by our eye classifi- 
cation plotted against that given in RC3 for 54 com- 
mon galaxies. The solid line connects points where 
our morphology matches with RC3 morpholgy. 

implying that our classifiers set the division criterion 
of E and SO to be more to the SO side than in other 
work ^\ Taking T < 1.5 as early type (E and SO) 
galaxies, the frequency of the early type galaxies in 
our sample 0.34 = 154/456 is consistent with what 
has been known from the work in the blue bands. In 
the following analysis we mostly take E and SO as a 
single class. To conclude this section, we quote the 
frequency of morphological types: 

E + SO : Sa : Sb : Sc : Sdm + Im = 0.34 : 0.18 : 0.21 : 0.21 : 0.06. 

(4) 

We do not detect any significant systematic change 
of morphological compositions with magnitude bins. 
So we regard this fraction to represent a fair value for 
field galaxies brighter than r* — 18 mag. 



In a later work, we tried to reclassify the sample by resetting 
the division criterion in a way that the frequency of E:S0 agrees 
with more traditional values in the blue band, but we do not 
adopt this reclassification in the present paper. 
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Fig. 4. — Distribution of morphological types for the 
456 galaxies (solid line). The dotted line and dashed 
line indicate the type mixes given in Loveday (1996) 
and Fukugita et al. (1998), respectively. 

4. Color indices 

4.1. Color indices as a function of morphol- 
ogy 

In Figure 5 we plot (a) u* - g*, (b) g* —r*, (c) 
r*—i*, and (d) i* — z* colors against T. The size ofthe 
symbol is again set so that its area is proportional to 
the number of galaxies falling in the grid. The mean 
and rms of the color distributions are calculated in Ta- 
ble 1, and represented in Figure 6 (a)-(d)0 As noted 
in section 2, colors of each galaxy are measured for a 
common aperture which is defined as 2 x rp(r*). The 
aperture is reasonably large, so that our color indices 
are regarded as representing nearly those for the to- 
tal flux. We correct all colors for Galactic reddening. 
The open squares connected by lines in Figures 5 and 
6 are colors calculated convolving Kennicutt's (1992) 

When calculating the mean, we dropped galaxies having too 
unrealistic colors; their colors are likely to be ascribed to errors 
due to unsuccessful deblending processes. 



Fig. 5. — Colors of our galaxies plotted against mor- 
phology. Panels (a), (b), (c), and (d) are for u' — g' , 
g' — r', r' — i', and i' — z', respectively. All colors 
have been corrected for Galactic reddening. The area 
of each filled circle is proportional to the number of 
galaxies falling into the grid. The open squares con- 
nected with solid lines indicate colors calculated (FSI) 
from Kennicutt's spectrophotomctric atlas of nearby 
galaxies. 

spectrophotometric atlas of nearby galaxies with the 
responses of the SDSS photometric system to give 
broad band colors (Fukugita, Shimasaku, & Ichikawa 
1995; hereafter FSI). 

We indicate in Fig. 6(a) representative uncertain- 
ties in our T measurements for early and late types 
by horizontal bars (the estimate at T = and 3) 
which are estimated from the difference between our 
classification and that of RC3 shown in Fig. 3. 

We observe (more clearly in Figure 6) that the data 
are distributed fairly close to the FSI calculations for 
(a) u* - g* , (b) g* — r*, and (c) r* - i* colors. For 



The uncertainty we plotted here is the quadratic sum of the 
random and systematic errors around the RC3 type at a given 
T. 
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Fig. 6. — Mean and rms colors of dcrcddcned galax- 
ies plotted against morphology. Panels (a), (b), (c), 
and (d) are for u' — g', g' — r', r' — i', and i' — z', 
respectively. The open squares connected with solid 
lines indicate synthetic calculation of FSI. The two 
filled squares with error bars in (a) show the errors in 
T(ours) measurements at T = and 3 estimated from 
the difference between our classification and that of 
RC3. 



g* — r* the difference between the mean of the data 
and the calculation is less than 0.05 except for Im (a 
0.3 mag discrepancy), for which the synthesis calcu- 
lation uses the bluest known galaxies (NGC4449 and 
NGC 4485) to define the envelope of the bluest color. 
The next largest discrepancy in g* — r* color is seen 
for elliptical galaxies: the data are redder than the 
synthesis calculation by 0.05 mag. This level of dis- 
crepancy, however, is seen for B — V color of elliptical 
galaxies, either with Kennicutt's spectrophotometric 
data or with the spectroscopic energy distribution 
(SED) of Coleman, Wu, & Weedman (1980, hereafter 
CWW), when compared with the broad band color: 
the template SED of elliptical galaxies is not suffi- 
ciently red to give B ~ V ~ 1 obtained by broad band 
photometry, as noted by FSI. 



For r* — i* color the difference between the data 
and the FSI synthetic colors is less than 0.04, again 
except for Im galaxies. The disagreement between 
our data and the FSI calculations is somewhat larger 
for u* — g* , but the scatter is wider for this color; we 
cannot conclude that the disagreement is significant. 
We remark, however, that the calculation being red- 
der than observation by 0.2 mag for E galaxies may be 
ascribed to the use of metal rich giant elliptical galax- 
ies in constructing the composite SED. The u* — g* 
index is most sensitive to the metal abundance. For 
Im galaxies the observed u* — g* color is redder than 
the calculation because of the use of a blue template 
object, as we have seen for the Im type in g* — r* . 

A significant deviation is seen for the i* — z* color, 
for which the data are clearly bluer than the syn- 
thetic calculation by 0.1 mag for early- type galaxies 
and 0.2 for later-type galaxies. While the SED that 
covers the z' band used in FSI relies on an extrapola- 
tion from the bluer wavelengths, we cannot blame this 
disagreement simply on an incorrect extrapolation, as 
we discuss in section 4.2. 

We present in Figure 7 color-color diagrams (a) 
g* - r* vs. u* - g* , (b) r* - i* vs. g* — r* , and (c) 
i* — z* vs. r* — i* . The dots represent each galaxy, 
and the curves are predictions from synthetic calcula- 
tions of FSI (filled squares), CWW SEDs extended to 
the near infrared by Neugebauer (personal communi- 
cation) (filled triangles), both of which use empirical 
SED of galaxies, Kodama & Arimoto's (1997) stel- 
lar population synthesis model (open squares), and 
GISSEL of Bruzual & Chariot's (1996, unpublished; 
see also idem 1993) stellar population synthesis model 
(open triangles). The convergence of the four syn- 
thesis calculations is very good for the g* — r* vs. 
u* - g* plot. A trend is visible that the SDSS pho- 
tometry gives a somewhat bluer u* — g* (or slightly 
redder g* — r*) for early type galaxies than all syn- 
thetic calculations. The agreement among synthetic 
calculations is still fairly good for the r*—i* vs. g*—r* 
plot; the SDSS data also closely follow the synthetic 
calculations. In the i* — z* vs. r* — i* plot, however, 
the data are clearly shifted by 0.1—0.2 mag down- 
ward (bluer i* — z*) compared with all four synthetic 
curves, which mutually agree very well. |^| The i' — z' 
color of the SDSS galaxies is deviated from all model 
calculations. 



20 This does not necessarily mean that model SED's are correct, 
since they are not well constrained beyond i' pass band. 
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Fig. 7. — (a) Color-color distribution of our galaxies: 
u* — g* vs g* — r* . The dots represent each galaxy, 
and the curves are predictions from synthetic calcula- 
tions of FSI (filled squares), CWW SEDs extended to 
the near infrared by Neugebauer (personal communi- 
cation) (filled triangles), Kodama & Arimoto's (1997) 
stellar population synthesis model (open square), and 
GISSEL of Bruzual & Chariot's (1993) stellar popu- 
lation synthesis model (open triangles). 

4.2. Comparison with other multicolor broad- 
band photometry 

Buta & Williams (1995, hereafter BW) present 
VRcIc colors for 501 nearby galaxies classified into 
morphology according to the RC3 system. UBV 
colors for these galaxies are given in RC3. In this 
subsection, we compare SDSS photometry with the 
BW (and RC3) catalog to examine the consistency 
among broad band photometries. Color transforma- 
tions from the U BVRcIc system to the SDSS system 
are carried out using the transformation laws given in 
FSI. The RC3 morphology indices are also converted 
into our system. We adopt BWs total color to com- 
pare with our Petrosian color. Figure 8(a)-(d) gives 
plots for the BW sample, similar to the plots in Figure 



Fig. 7. — (b) Same as (a), but for g* - r* vs r* — i* . 

5(a)-(d). We remark that the BW sample is not ho- 
mogeneous with respect to morphology, but is highly 
biased towards early type galaxies. We expect, how- 
ever, that this does not produce any biases to our 
color analysis. We overlay the data of our sample 
(Fig. 5; we omit T = 0.5 and 1.5 data since no galaxy 
is classified as these types in the BW sample) by open 
circles, with a horizontal offset of +0.2 for clarity. 

A comparison of Figure 8 with 5 shows that the 
u' — g', g' — r' , and r' — i' distributions of the BW 
galaxies are grossly consistent with those for SDSS 
photometry. The agreement of the u' — g' colors be- 
tween the two samples is particularly good, although 
a somewhat smaller scatter for the early type galax- 
ies of the BW sample may be ascribed to its selection 
effect that the BW sample is not homogeneous. The 
g' — r' and r' — i' colors for early type galaxies are 
slightly bluer for the BW sample. The appreciable 
discrepancy found in the comparison of the two sam- 
ples is with i' — z' colors, i' — z' colors of the SDSS 
sample are systematically bluer than those of BW by 
~ 0.1 — 0.2 mag. This trend is similar to what was 
found between the FSI calculation and the SDSS mea- 
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Fig. 7. — (c) Same as (a), but for r* — i* vs i* — z*. 

sures of galaxies in the previous subsection. 

This discrepancy in i' — z' color leads us to sus- 
pect that there may be a systematic error by 0.1 — 0.2 
mag in z' band photometry, but the real reason is 
not clear to us. The preliminary work for the SDSS 
primary standard stars (Smith ct al. 2001) shows 
that all SDSS color indices relative to the Johnson- 
Morgan-Cousins photometric system differ from a 
synthetic calculation using Gunn-Stryker stars (Gunn 
& Stryker 1983) by no more than 0.03 mag. We 
have recently recognized some unexpected shifts of 
the filter response function in the 2.5 m camera sys- 
tem due to the evaporation of water vapor from sil- 
icon layer of the coating surface when the filters are 
soaked in a vacuum chamber, but we estimate the 
effect to be smaller than 0.015 mag for all color 
bands. The difference in the response between the 
front-illuminated thick CCD device used in the 2.5m 
camera and the back-illuminated thin device for the 
Photometric Telescope (which defines the SDSS mag- 
nitude) does not cause an offset of more than 0.01 
magnitude. The sky brightness and its time varia- 
tion are the highest in the z' band, so errors in the 



Fig. 8. — Same as Fig. 5, but for Buta & Williams' 
(1995) galaxies. Overlaid are the plots for our sample 
(Fig. 5) by open circles with a horizontal offset of +0.2 
for clarity. 

sky subtraction in z' band could affect photometry of 
galaxies. However, it is not easy to estimate photo- 
metric errors due to the sky subtraction errors, since 
they affect galaxy photometry in a complicated way 
depending on sizes and surface brightnesses of galax- 
ies. At present, we have not found clear evidence that 
they account for the discrepancy in i' — z' color. It 
should be, however, worth noting that the observed 
amount of the discrepancy in i' — z' does not depend 
either on the brightness (and size) of galaxies nor on 
the position of galaxies on a strip, ie, when they were 
imaged. The effect of 0.1-0.2 mag seems a bit too 
large to be ascribed to the uncertainties in the present 
photometric system, but we cannot find the reason for 
this discrepancy at the moment. 

5. Scale lengths 

In this section, we study the relation between ap- 
parent size and apparent brightness of galaxies. We 
expect different relations for different morphologi- 
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cal types. The relation may also depend on the 
color band used, if star forming regions or metallic- 
ity distributions are not homogeneous throughout the 
galaxy. The information on the radius-brightness re- 
lation is also needed to construct a realistic simulation 
of galaxies. We additionally examine whether the re- 
lation derived from SDSS photometry is consistent 
with our knowledge of the galaxy sizes we obtained 
from other studies of local galaxies. The consistency 
would allow us to supplement the SDSS data with 
the results from other, more detailed studies, such as 
which have not yet been done with SDSS. 

We consider r 50 estimated in the five color bands 
as a measure for the galaxy size. Note that r 50 is 
approximately the conventionally defined effective ra- 
dius (half light radius). We adopt Petrosian magni- 
tudes, mp, also in the five color bands to represent 
apparent brightness. 

Four panels in Figure 9 show the plot of rso(Aj) 
(Aj = u' , g' , r' , and i') as a function of apparent 
brightness in the r' passbands. We divide our galaxy 
sample into three morphological classes on the basis 
of T values: the open circles are early type galaxies 
with T < 1.5 (E and SO galaxies), the filled circles 
denote 1.5 < T < 4 (SO/a to Sbc), and the crosses 
are for late spiral galaxies with T > 4 (later than Sc). 
The data are well represented by the curve 

logr 50 (Aj)" = -0.2 [mp(Xj) - 16 mag] + ay, (5) 

where ay are constants that depend on Aj and Xj. 
This relation holds for any combination of the pass- 
bands Xi,Xj with the parameters ay (and their rms 
scatters) being obtained by least squares fits to the 
plots for the three morphological classes. These pa- 
rameters are given in Table 2 for Aj = u' , g', r' , i' , z', 
Xj = r', and for Aj = Aj = u' , g' , r', i', z'. 

The plots in Figure 9 show that half light Pet- 
rosian radii r§o are nearly independent of the color 
bands with which these radii are calculated. Looking 
into details, however, we see the trend that the half 
light radius shrinks slightly from UV to redder bands 
for all morphological types (see rows 1-5 of Table 2). 
For early type galaxies we observe a 10% decrease in 
the effective radius continuously from the v! to the 
z 1 band. This is interpreted as due to the effect of 
the color gradient of early type galaxies driven from 
the metallicity gradient. For spiral galaxies the de- 
crease is 20-25%. This represents the star forming 
component being more extended towards outer parts 
of galaxies: in other words an increasingly dominant 
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Fig. 9. — (a) Half light Petrosian radii measured in u', 
plotted as a function of Petrosian magnitude in the 
r' band. Magnitudes have been corrected for Galac- 
tic extinction. Open circles, filled circles, and crosses 
indicate T < 1.5, 1.5 < T < 4, and T > 4, respec- 
tively. The solid, dotted, and dashed lines in each 
panel show the best fit of eq. (5) to the three mor- 
phological classes of galaxies, respectively. 

importance of the bulge component in redder pass- 
bands. 

The fitting formulae with i = j give useful param- 
eters (see rows 6-10 of Table 2) to construct simu- 
lations of galaxies. The knowledge of these param- 
eters is particularly important in understanding the 
detection limit of galaxies and the performance of star 
galaxy classification. 

Let us now investigate the consistency of eq.(5) 
with the relation involving the effective radius known 
in the literature. The most accurate study concern- 
ing the scale of galaxies is presented by Kent (1985) 
on the basis of excellent Thuan-Gunn r band CCD 
images of nearby bright galaxies with known red- 
shifts. He derived the relations between intrinsic size 
and surface brightness separately for the bulge and 
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Fig. 9. — (b) Same as (a), but for g' . 

disk components. Using the transformation law of 
FSI, we convert his relations between the effective ra- 
dius (kpc) and the effective surface brightness (mag 
arcsec -2 ) to those in the SDSS r' band as, 

[i e = 2.5 log[r e (kpc)] + 20.3(mag sec" 2 ) , (E, SO) 

/i e = 2.51og[r e (kpc)] + 20.5(mag sec 2 ) , (disks) 

(7) 

where we have assumed the Hubble constant to be 
100 km s _1 Mpc" 1 . 

When these relations are fitted to the SDSS data 
for T < 1.5 (bulge dominated) and T > 4 (disk 
dominated), we obtain for the effective radius r e de- 
fined by J(r) oc exp[-0.7289(r/r e ) 1 / 4 ] for the de 
Vaucouleurs profile and the scale length defined by 
7(r) oc exp(— r/h) for the exponential disk, 

r e = 2.561q g2 kpc , (de Vaucouleurs profile) (8) 
and 

h = 2.62±°;g3 kpc , (exponential disk) (9) 




11 12 13 14 15 16 17 18 
m p (r*) (mag) 



Fig. 9. — (c) Same as (a), but for r'. 

where the errors indicate the la range. Note that 
the average luminosity of galaxies that contribute to 
galaxy number counts in the bright magnitude is very 
close to L* (within ±0.1 mag error for r* = 15 — 18 
mag) with the slope of the Schechter function a = 
— 1.2. Therefore, the values given in eqs. (||) and (||) 
represent the average sizes for the L* galaxies. 

Our data indicate that the mean h derived from 
late type disks is about the same value as the mean 
r e of early type galaxies. With these two parameters, 
our data are consistent with the relations derived by 
Kent (1985). 

6. Concentration index 

We define the (inverse) concentration index by the 
ratio of the two Petrosian radii C = r$o /rgo measured 
in the r' band. Figure 10 shows C against T (see 
also Table 3). We have dropped 30 galaxies having 
r 5Q < 2", since large smearing effects are suspected 
for C for such small galaxies unless significant cor- 
rections are made for the seeing (average FWHM is 
~ l."5). The two filled squares with error bars near 
the top of the panel show the typical errors in T(ours) 
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Fig. 9. — (d) Same as (a), but for i'. 

measurements at T — and 3. 

A tight correlation is seen between C and T; C in- 
creases with T, which is understood by the spheroid 
or bulge dominance of early type galaxies and an in- 
creasing disk dominance of later- type galaxies. We 
have studied similar plots for galaxies with magni- 
tudes divided into three bins, r* < 16, 16 < r* < 17, 
and 17 < r* < 18. No systematic effects are observed 
among the samples with different apparent bright- 
nesses (hence different apparent sizes) of galaxies, in 
so far as we limit ourselves to galaxies with rso > 2". 
The galaxies with r^o < 2" show a weaker correlation 
between C and T due to smearing by seeing but also 
possibly due to errors of eye classification. [^] 

The scatter in C for a given T is larger than ex- 
pected from the errors in our T measurements. This 
reflects the fact that galaxies with a given morphol- 

We have examined whether highly inclined galaxies would dis- 
turb the correlations between C and T. We have found that the 
correlation does not change appreciably againt the sample with 
specified minor-to-major axis ratio. Our concentration index 
defined using the Petrosian radii is a robust indicator against 
inclination of galaxies. 



Fig. 10. — Concentration index C plotted against 
eye morphology T for 426 galaxies having r$o > 2". 
The vertical line indicates the boundary which divides 
E/SO galaxies and spiral/Im galaxies on the basis of 
eye morphology. The horizontal line is for the bound- 
ary accepted in the text {C\ — 0.33) dividing low-C 
('early type') and high-C ('late type') galaxies. The 
two filled squares with error bars near the top of the 
panel show the errors in T(ours) measurements at 
T = and 3 estimated from the difference between 
our classification and that of RC3. 

ogy have a fairly wide range of the bulge-to-disk ratio 
(e.g., Kent 1985; Simien & de Vaucouleurs 1986; Ko- 
daira, Watanabe, & Okamura 1986) and the bulge- 
to-disk ratio is directly related to the concentration 
index. 

We examined the correlations of visual morphol- 
ogy with a number of parameters measured by Photo 
(likelihoods of a fitting to the de Vaucouleurs or expo- 
nential profiles, texture parameters defined by Photo) 
in addition to the three parameters studied in this pa- 
per. The correlation of the likelihood of profile fitting 
and of the texture parameter with our visual morphol- 
ogy is not strong (we used the version of Photo version 
5_0_3 of late 1999). The concentration index shows 
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the strongest correlation with visual morphology, and 
a combination of the concentration index with other 
parameters, such as surface brightness, color, asym- 
metry parameter, does not appreciably enhance the 
correlation. 

The concentration index is perhaps the best pa- 
rameter to be used to classify morphology of galax- 
ies. This is basically the same conclusion as reached 
by Doi et al. (1993), and further corroborated by 
Abraham et al. (1994). These authors used a com- 
bination of a concentration index and mean surface 
brightness, while the concentration index is sufficient 
for our case. This simplification arises merely from 
the different definition of the concentration indices. 
In Doi et al. and Abrahams et al. the concentration 
index is defined using isophotal photometry. The cal- 
culated index then depends on both apparent bright- 
ness and distance of galaxies. To correct for these re- 
dundant effects the second parameter is needed, and 
those authors took average surface brightness, which 
also depends on apparent brightness and the distance. 
In our case the Petrosian quantities are physically well 
defined, independent of the apparent brightness or the 
distance. Therefore the introduction of a second pa- 
rameter is unnecessary. 

The correlation in Figure 10 suggests that the C 
parameter can be used for morphological classifica- 
tion allowing for the uncertainty of roughly AT ~ 1.5 
(see Table 3 above). The most important problem 
in the classification using the C parameter is that 
the C parameter for Sa galaxies can be as small as 
that for E galaxies, making it difficult to construct a 
purely early type galaxy sample free from a contam- 
ination of Sa galaxies. As noted above this problem 
is not ameliorated with the use of any second param- 
eters we studied. On the other hand, the distinction 
between SO and Sa is not too difficult, allowing for 
some intermediate cases, upon visual inspections due 
to the presence of spiral arms and/or HII regions in 
Sa galaxies, which are not properly picked up with 
the parameters given by Photo. 

A practically important application is to distin- 
guish E/SO galaxies from spiral galaxies. Here we 
consider this problem in a more quantitative way, 
and study the completeness and contamination of the 
morphologically classified sample with the use of the 
C parameter. The upper panel of Figure 11 shows the 
completeness as a function of C: the solid curve rep- 
resents the completeness of early type galaxies when 
sample is selected with C < C\ for a given Ci, while 




0.2 0.25 0.3 0.35 0.4 0.45 0.5 



Fig. 11. — Completeness and contamination of early 
and late type galaxies as a function of C\ . The solid 
and dotted lines are for early and late type galaxies, 
respectively. 

the dashed curve represents the completeness of the 
late type galaxy sample with C >C\. The complete- 
ness of the two samples balances with C\ — 0.35 at 
80%. On the other hand, the lower panel of Figure 11 
indicates the contamination from the opposite type. 
Namely, the solid curve is the contamination from 
late type galaxies to the early type galaxy sample 
with C < C\. The dashed curve is the contamination 
by early type galaxies to the late type sample with 
C > C\. The two curves become even at C\ = 0.32. 

The inspection of two curves in the lower panel of 
Figure 11 shows that we can make a reasonably pure 
late type galaxy sample with the choice of large C\ , 
say, 0.37 at a loss of completeness (70% complete- 
ness; contamination is less than 5%). The opposite 
is not true: one cannot decrease the contamination 
of an early type sample by late type galaxies below 
15%, whatever value of C\ is chosen, as a result of a 
contamination from Sa galaxiespl 

We suspect that this contamination is not due to Seyfert galax- 
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Our recommended choice of C\ is 0.33 if we want 
dichotomous classification of galaxies into early and 
late types. For this specific choice the completeness 
and contamination are given in Table 4. 93 out of the 
136 early type (E/SO) galaxies are classified as 'early' 
types, i.e., a completeness of 68%. On the other hand, 
24 'early'-type galaxies are not real E/SO galaxies, 
which means a contamination of 21% (24/117). Simi- 
larly, the fraction of late type (Sa-Im) galaxies classi- 
fied as 'late' types is 92% (266/290) and the contam- 
ination in the 'late' galaxy sample from early type 
galaxies is 14% (43/309). 

In summary, our analysis demonstrates that the 
parameter C — r 50 / r 90 is useful for automated classi- 
fication of early and late type galaxies, if one is satis- 
fied with a completeness of ~ 70 — 90%, allowing for 
a contamination of ~ 15 — 20%. It appears difficult to 
improve the performance with the use of simple pho- 
tometric parameters. Further improvement requires 
the introduction of more complicated parameters such 
as those characterizing the presence of spiral arms 
and HII regions. An improvement of late type galaxy 
classification may be achieved with the introduction 
of the asymmetry parameter (e.g., Abraham ct al. 
1996). 

7. Conclusions 

We have produced a morphologically classified sam- 
ple of 456 bright galaxies from imaging data which 
were taken during the Sloan Digital Sky Survey com- 
missioning phase. We have classified galaxies into 
seven Hubble types (E, SO, Sa, Sb, Sc, Sdm, and Im) 
by visual inspections carried out by four individuals. 
We have measured the distributions of colors, scale 
lengths, and concentration indices in the SDSS photo- 
metric system as functions of morphology. Our prime 
purposes are (i) to derive the fundamental statisti- 
cal properties of galaxies in the SDSS passbands, (ii) 
to examine the consistency with existing galaxy data 
other than the SDSS, and (iii) to find the parameters 
that can be used for an efficient automated morpho- 
logical classifier. All parameters we have studied in 
this paper, except for visually determined morphol- 
ogy, are among the standard outputs from the SDSS 
image analysis software (Photo), and thus routinely 
registered in the SDSS photometric catalog. 

We have found that the colors of our galaxies in 

ies, since the color of the contaminants are not bluer than 
others. 



the SDSS system match with those obtained from 
template SED data and from UBVRcIc broad band 
photometry of nearby galaxies to within «0.05 mag. 
The exception is i' — z' color for which the discrep- 
ancy is as large as 0.1-0.2 mag. We could not find the 
origin of this discrepancy yet. 

We found that the half light radius falls on the line 
logr5o = —0.2m + const, and the constant we found 
indicates r e ~ 2.6 kpc for de Vaucouleurs profiles, 
and h ~ 2.6 kpc for exponential disks of late type 
spiral galaxies with the aid of the surface brightness 
radius relation derived by Kent (1985). We found 
additionally that the half light radius of galaxies de- 
pends slightly on the color band and the dependence 
is consistent with the expected distribution of star 
forming regions for late type galaxies and with the 
known color gradient in early type galaxies. 

Finally, we have shown that the concentration in- 
dex shows a strong correlation with the morphological 
types. We explored the possibility of using this pa- 
rameter for automated morphological classification of 
galaxies, in order to classify a huge number of galax- 
ies into early, intermediate, and late types with pa- 
rameters that are routinely produced in the photo- 
metric data processing of the SDSS. For an optimum 
ratio of [r 5 o/r 9 o]i = C\ = 0.33, we find, however, a 
15 — 20% contamination from opposite types of galax- 
ies. The contamination is particularly important for 
the early type galaxy sample, and we must consider 
some higher order photometric parameters to enhance 
the performance. 

Note Added 

Photo has been continuously updated over the 
years. After we submitted our paper, Strateva et al. 
(submitted to AJ) found on the basis of a different 
sample (287 galaxies at g* < 16) that the likelihood 
of profile fitting based on a later version of Photo (5_2) 
correlates with eye morphology better than we found. 
However, they also found that the profile fitting does 
not work well in comparison with concentration in- 
dex for bright galaxies such as those we study in our 
paper. 

We thank the referee, Roberto Abraham, for his 
valuable comments. 

The Sloan Digital Sky Survey (SDSS) is a joint 
project of The University of Chicago, Fermilab, the 
Institute for Advanced Study, the Japan Participa- 
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Planck-Institute for Astronomy (MPIA), the Max- 
Planck-Institute for Astrophysics (MPA), New Mex- 
ico State University, Princeton University, the United 
States Naval Observatory, and the University of Wash- 
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Table 1: Mean colors of SDSS galaxies as a function 



of T. 



T 


Hubble type 


u* - 


Q* 


Q* - 


r* 


r* - 


i* 


i* - 


z* 


T < 0.5 


E 


1.79 


(0.26) 


0.83 


(0.14) 


0.41 


(0.05) 


0.27 


(0.06) 






86 




87 




87 




87 




0.5 < T < 1.5 


SO 


1.66 


(0.30) 


0.75 


(0.14) 


0.38 


(0.05) 


0.26 


(0.06) 






67 




67 




67 




67 




1.5 < T < 2.5 


Sa 


1.49 


(0.32) 


0.68 


(0.15) 


0.38 


(0.07) 


0.25 


(0.08) 






81 




81 




81 




81 




2.5 < T < 3.5 


Sb 


1.40 


(0.28) 


0.62 


(0.13) 


0.35 


(0.08) 


0.20 


(0.09) 






95 




94 




93 




94 




3.5 < T < 4.5 


Sc 


1.28 


(0.33) 


0.46 


(0.12) 


0.27 


(0.09) 


0.07 


(0.16) 






94 




98 




97 




91 




4.5 < T < 5.5 


Sdm 


1.16 


(0.38) 


0.42 


(0.24) 


0.20 


(0.12) 


0.04 


(0.20) 






22 




22 




21 




18 




5.5 < T 


Im 


1.14 


(0.33) 


0.61 


(0.23) 


0.23 


(0.08) 


0.05 


(0.17) 






6 




6 




5 




5 





Note. — The numbers in the parentheses are dispersions. The numbers in the second rows for each morphology listing arc the number 
of galaxies used to calculate the mean. Only galaxies satisfying the following criteria are used to calculate the mean: < u* — g* < 2.5, 
< g* - r* < 1.5, < r* - V < 1, -0.5 < V - z* < 0.5. 



Table 2: Fit coefficients aij in the relation between 
apparent half light radius and apparent brightness. 



Xi X 3 All T<1.5 1.5<T<4 T>4 



1 


u' 


r 1 


0.51 


(0.19) 


0.37 


(0.15) 


0.56 


(0.16) 


0.63 


(0.20) 


2 


g' 


r' 


0.51 


(0.18) 


0.36 


(0.12) 


0.55 


(0.14) 


0.68 


(0.17) 


3 


r' 


r' 


0.49 


(0.18) 


0.36 


(0.12) 


0.52 


(0.14) 


0.66 


(0.17) 


4 


i' 


r' 


0.48 


(0.17) 


0.35 


(0.12) 


0.51 


(0.14) 


0.64 


(0.16) 


5 


z' 


r' 


0.44 


(0.16) 


0.33 


(0.12) 


0.47 


(0.14) 


0.57 


(0.16) 


6 


u' 


u' 


0.94 


(0.16) 


0.87 


(0.15) 


0.96 


(0.14) 


0.99 


(0.17) 


7 


g' 


g' 


0.64 


(0.16) 


0.52 


(0.13) 


0.67 


(0.13) 


0.77 


(0.16) 


8 


r' 


r' 


0.49 


(0.18) 


0.36 


(0.12) 


0.52 


(0.14) 


0.66 


(0.17) 


9 


i' 


i' 


0.41 


(0.18) 


0.27 


(0.12) 


0.44 


(0.14) 


0.59 


(0.17) 


10 


z' 


z' 


0.34 


(0.19) 


0.20 


(0.12) 


0.36 


(0.16) 


0.53 


(0.19) 



Note. — See eq. (5) for the definition. The numbers in parentheses are dispersion around the fit. 
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Tabic 3: Concentration index C = r 5a /r 90 . 



T(our) 


N 


C 

mean 


rms 


T < 0.5 


80 


0.317 


0.027 


0.5 < T < 1.5 


56 


0.331 


0.037 


1.5 < T < 2.5 


73 


0.362 


0.054 


2.5 < T < 3.5 


92 


0.404 


0.057 


3.5 < T < 4.5 


97 


0.467 


0.052 


4.5 < T < 5.5 


22 


0.452 


0.055 


5.5 < T 


6 


0.439 


0.061 



Table 4: Correlation between eye morphology and 
that classified with the concentration index with C\ = 
0.33. 





T < 1.5 


T > 1.5 


sum 


contamination 


'early type' galaxies 


93 


24 


117 


21 % 


'late type' galaxies 


43 


266 


309 


14% 


sum 


136 


290 


426 




completeness 


68% 


92% 
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